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BACKGROUND OF THE INVENTION 



1 . Field of the Invention 

The present invention relates to processor synchronization, 
and more particularly to clock- synchronization in distributed 



3 2 • Discussion of Related Art 

J l j Quartz crystals are known to provide accurate and constant 

clocks for electronic equipment. For this reason they are widely 

used in relatively inexpensive electronic equipment such as PCs, 

fyL5 mobile phones, wrist watches, etc. Typically, quartz crystals 

80 are grown synthetically. The frequency of a quartz crystal is 

N 1 " determined by its thickness, d, according to: 



10 systems. 



f[kHz] = 



N[kHz*mm] 



Equation 1 



d[mm] 



20 



For the most common AT-cut, the constant N is N=1660 



kHz*mm. AT specifies a specific plane relative to the crystal 



axes of the quartz. The frequency ranges from 800khz up to 360 



Mhz, and the fundamental frequency ranges up to about 40Mhz. 



25 



The advantage of the AT-cut in comparison to other cuts is 



that the resonance frequency of the AT-cut is substantially 
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independent of temperature. It follows a third grade equation 
(see Fig. 1, and Equation 2) with a turning point, according to 
the form of the quartz and frequency at 25 to 30 degrees 
Celsius . 



Af T \,„( T T y Equation 2 



with a 3 =1.05 -10' 



a, = 0.0085 • L<p 



yjLO A^ = ^ 2Z -^? 0 in angle minutes, 

s in ppm and T t T inv in degrees Celsius. 



The gradient in the turning point is determined by the cut 
angle <p zz . A<p is the difference between the zero angle (p 0 at 
15 which the curve has a horizontal tangent in its turning point 
(see Fig. 1) . 

There are two main sources of frequency inaccuracies for 
commercial quartz oscillators: manufacturing inaccuracies and 
thermal effects. Manufacturing inaccuracies result from the 
2 0 finite mechanical accuracy of quartz thickness, as can be seen 
from Equation 1, and accuracy of the cut angle (see Equation 2) , 
which can only be adjusted to 10 angle minutes. Inexpensive 
quartz oscillators have a typical relative accuracy better than 



+/-10" 4 , which is adequate in most cases. As an example, most ITU 
modem recommendations specify a clock accuracy of +/-10" 4 . In 
inexpensive equipment, this is the main source of frequency 
offset . 

5 As shown in Fig. 1, additionally the quartz frequency 

depends on the temperature of the environment. In a room 
temperature environment of 15°C ... 25°C, the thermal inaccuracy 
M* of quartz oscillators is insignificant compared to manufacturing 

JKSSi 

O inaccuracy. Further, there are still other reasons for frequency 
Wlo inaccuracies of quartz oscillators such as long term drift and 
S = J phase noise. However, these inaccuracies are typically 

negligible compared to the above effects. 
y : There are many ways to improve the accuracy of oscillators, 

O such as parallel adjustable capacitors or mounting in a 
15 temperature regulated housing. However, oscillators in standard 

PCs do not use such sophisticated techniques to limit complexity 

and costs. 

One area where clock inaccuracies can be particularly 
troublesome is in distributed computing. The clocks of different 
2 0 processors need to be synchronized to limit errors. 

Synchronization can be a particularly difficult problem. For 
example, if the clocks of collaborative video discussion par- 
ticipants' are to be synchronized to a desirable level, e.g., up 
to 3 0 ms, it is not enough to adjust local clocks only at the 
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beginning of a collaborative video discussion session. The 
clocks need to be monitored and adjusted continuously. 
Otherwise, the clocks will drift apart because of, for example, 
temperature drift, aging and the mechanical accuracy at which 
5 the frequency can be set by cutting. The drift rate is the 

deviation from the Coordinated Universal Time (UTC) , divided by 
the time T, the measuring period. 
N : One approach to distributed clock synchronization was 

reported by Flaviu in 1989 in his article "Probabilistic Clock 
SLO Synchronization' 7 , which is the basis for many synchronization 
SI algorithms. It is called probabilistic, because in the approach 
L it is not guaranteed that a client can always read the clock of 
U a server with an a-priori defined precision. By means of several 
O successive reading attempts, the client is able to synchronize 
15 onto the server's clock with an adjustable precision, and with a 
probability of success which is also adjustable via a maximum 
number of allowed reading attempts. The disadvantages of this 
method include: the method does not take any processing time 
into account, and no dynamic adjustment to the current traffic 
20 situation in a TCP/IP network is implemented. Another 

disadvantage of the method include to achieve a good precision, 
a large load may be placed on the network, caused by the syn- 
chronization messages. 
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Another clock synchronization method is used in the Network 
Time Protocol (NTP) , which is an Internet Standard Recommended 
Protocol, described in the RFC-1119, RFC-1305 and RFC-2030. For 
the basic approach to measure the offset of another computer 
over TCP/IP, NTP is uses substantially the same algorithm as 
presented in David L. Mills, RFC-2 03 0 Simple Network Time 
Protocol (SNTP) Version 4 for IPv4, IPv6 and OS I, Network 
Working Group, October 1996. To save messages, and therefore 
reduce the network load caused by the update messages, an 
additional "logical clock" is implemented on the client side. 
This logical clock controls the local hardware clock, so that 
high precision can be realized, without many synchronization 
meS sages. Although this protocol works well, and high precision 
can be normally achieved (up to a few tens of milliseconds in 
global WANs), a major disadvantage of NTP is its slow speed. 
Based on the experience of the Rutgers University with NTP, it 
may take up to 30 minutes for the client's clock to synchronize 
for the first time on a time server. If, in addition, the 
difference between the two clocks is more than a few minutes, it 
can take much more time for synchronization to occur for the 
first time. 

Other minor disadvantages of NTP include: the overhead of 
the NTP-messages, because of security and redundancy features, 
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and the synchronization is typically done on more than one 
server, which is until now was neither possible nor necessary. 

Therefore, a need exists for a system and method for 
increasing the precision of clock synchronization in distributed 
5 systems . 

SUMMARY OF THE INVENTION 

jf According to an embodiment of the present invention, a 

method is provided for synchronizing distributed processors. The 
^LO method comprises determining a desired number of offset values 
nt between two processors, wherein each processor comprises a 

quartz crystal, determining parameters of a regression line, 
wherein the regression line is a function of the offset values 
0 over the desired number of offsets, and adjusting a 
15 synchronization interval according to the parameters. 

The desired number of offset values is thirty. 
The offset values are functions of a difference between 
relative inaccuracies corresponding to the quartz crystal of 
each processor. 

20 Determining the parameters further comprises fitting a 

straight line, y = a + jb • x to a collection of N measurement 
pairs [yi/Xi) with minimum mean square error, wherein a and b are 
the parameters. 
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According to an embodiment of the present invention, a 
program storage device is provided readable by machine, tangibly 
embodying a program of instructions executable by the machine to 
perform method steps for synchronizing distributed processors. 
The method comprises determining a desired number of offset 
values between two processors, wherein each processor comprises 
a quartz crystal, determining parameters of a regression line, 
wherein the regression line is a function of the offset values 
over the desired number of offsets, and adjusting a 
synchronization interval according to the parameters. 

According to another embodiment of the present invention, a 
method is provided for synchronizing distributed processors. The 
method comprises establishing a socket -connection between at 
least two processors, determining a roundtrip delay, and 
determining a roundtrip-delay threshold. The method further 
comprises determining a current round- trip delay and an offset, 
adding the current round-trip delay to a list of roundtrip 
delays, and determining a new roundtrip-delay threshold. The 
method comprises determining whether the current roundtrip delay 
is greater than the new threshold: upon determining the current 
roundtrip delay to be greater than the new threshold, 
determining whether a desired number of round-trip delays have 
been determined; and upon determining that the current threshold 
is not greater than the new threshold, determining whether the 



offset is greater than an offset threshold. The method further 
comprises adjusting a clock according to an offset, and 
determining a linear regression. 

A probability of the round-trip delay being greater than 
the roundtrip-delay threshold is about 0.5 and a probability of 
the round-trip delay being less than the roundtrip-delay 
threshold is about 0.5 

Determining whether thirty round- trip delays have been 
determined further comprises entering a synchronization method 
upon determining the desired number round-trip delays. 
Determining whether thirty round-trip delays have been 
determined further comprises determining a current round- trip 
delay and an offset upon determining less than the desired 
number delays. 

Adjusting a clock according to an offset further comprises 
decrementing by an update- interval upon determining the offset 
to be greater than the offset threshold, and incrementing by the 
update- interval upon determining the offset to be less than the 
offset threshold. 

The method comprises determining, recursively, a current 
round-trip delay and an offset. 

Determining a linear regression further comprises setting a 
current synchronization time. The method comprises determining 
whether a number of measured offsets is greater than a desired 



number: upon determining that the number of offsets is greater 
than the desired number, removing an oldest offset from a list 
of offsets and adding a current offset to the list and 
determining parameters of a regression line from the list of 
5 offsets, and upon determining that the number of measured 
offsets is not greater than the desired number, adding the 
current offset to the list. The method further comprises 
estimating the current offset using the regression line, and 
JrJ incrementing the current synchronization time. The method 

comprises determining whether the current synchronization time 
m is greater than an update- interval : upon determining the current 
U synchronization time to be less than the update - interval , 
hi estimating the current offset using the regression line, and 
O upon determining the current synchronization time to be greater 
15 than the update - interval , measuring a current roundtrip delay 
and offset. 

The desired number of roundtrip delays is thirty. 

According to an embodiment of the present invention, a 
system is provided for synchronizing distributed processors. The 
2 0 system comprises a first processor connected to a network, 
wherein the first processor sends a sync-request message 
comprising a time current local time of the first processor. The 
system further comprises a second processor connected to the 
network and connected to the first processor via the network, 
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wherein the server receives the sync -request message, and stores 
a time of arrival of the sync-request message and sends a sync- 
response message the first processor, wherein the sync-response 
message comprises the current local time of the first processor, 
5 the time of arrival and a current local time to the second 
processor. 

M» BRIEF DESCRIPTION OF THE DRAWINGS 

M Preferred embodiments of the present invention will be 

JjjLO described below in more detail, with reference to the 
K accompanying drawings: 

L Fig. 1 shows a typical drift for various AT crystal cut 

M angels; 

□ Fig. 2 is a distribution of roundtrip delays; 

15 Fig. 3 is a roughly outlined histogram of the medium 

period; 

Fig. 4 is a roughly outlined histogram of the rare period; 
Fig. 5 is a roughly outlined histogram of the busy period; 
Fig. 6 shows messages for synchronization over TCP/IP, 
2 0 according to an embodiment of the present invention ; 

Fig. 7 shows the total offset between two unsynchronized 
computers, according to an embodiment of the present invention; 
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Fig. 8 is a transient response of the interval between two 
synchronization attempts using linear regression, 19.6kbps 
connection, according to an embodiment of the present invention; 

Fig. 9 illustrates an algorithm with linear regression, 
client side, according to an embodiment of the present 
invention; 

Fig. 10 is an algorithm of the sub function l_regression, 
according to an embodiment of the present invention; 

Fig. 11 is a transient response of the interval between two 
synchronization attempts using linear regression, 19.6kbps 
connection, according to an embodiment of the present invention; 
and 

Fig. 12 shows a measured offset between two synchronized 
computers (using linear regression), connected via a 19.6 kbps 
modem, according to an embodiment of the present invention. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

It is to be understood that the present invention may be 
implemented in various forms of hardware, software, firmware, 
special purpose processors, or a combination thereof. In one 
embodiment, the present invention may be implemented in software 
as an application program tangibly embodied on a program storage 
device. The application program may be uploaded to, and executed 
by, a machine comprising any suitable architecture. Preferably, 
the machine is implemented on a computer platform having 
hardware such as one or more central processing units (CPU) , a 
random access memory (RAM), and input/output (I/O) interface (s) . 
The computer platform also includes an operating system and 
micro instruction code. The various processes and functions 
described herein may either be part of the micro instruction 
code or part of the application program (or a combination 
thereof) which is executed via the operating system. In 
addition, various other peripheral devices may be connected to 
the computer platform such as an additional data storage device 
and a printing device. 

It is to be further understood that, because some of the 
constituent system components and method steps depicted in the 
accompanying figures may be implemented in software, the actual 
connections between the system components (or the process steps) 
may differ depending upon the manner in which the present 
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invention is programmed. Given the teachings of the present 
invention provided herein, one of ordinary skill in the related 
art will be able to contemplate these and similar 
implementations or configurations of the present invention. 
5 In a scenario where the content of a video needs to be 

discussed by two or more participants in a teleconference 
session in real time, it is important that video is delivered 
H= synchronously to all participants. Without collaborative video 
S delivery, it may be difficult to discuss the content, as the 
^LO content would be different at every participant. 
Jf| One solution would be to synchronized the clocks up to at 

L least 30 ms, similar to the number of frames per second in the 
U> MPEG standard, namely 2 9.97 frames per second. Thus, the 
O duration of each frame can be computed to 33.37 ms, which is 
15 roughly the 30 ms mentioned above. If the clocks are not 

synchronized, the quality of the synchronized video presentation 
can deteriorate. For example, if a camera recording video is 
moved quickly, or there is a sudden switch from one source 
channel to another. If in this moment, a client pauses the video 
20 stream and the clocks are not synchronized, each client will see 
at least one more frame. The server can then send the current 
frame to each client, but it can take an undesirable amount of 
time until the frames reach the clients. 
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By synchronizing a starting point of the video playback 
between all participants and the server to be an absolute time 
in which the server tells the clients when to start the 
playback, all the timestamps of the following frames are 
referenced to the updated starting point. A relative timestamp 
such as, "Start playing frame number x in 30 seconds" may not be 
useful because of unknown network delays, which cannot be 
estimated on 30 ms accuracy. Network delay is a measurement of 
the time it takes a message traveling from the server to reach 
the client. Further, local internal system clocks on local 
devices may have been set manually to within a minute or two of 
the actual time and are rarely reset at regular intervals. 
Finally, the present invention can take into account time 
differences between different time zones. 

In distributed systems the task of synchronizing clocks is 
made difficult, among other reasons, because of the existence of 
hardly predictable communication delays. There is an arbitrary, 
random time delay between a time when a message is sent from one 
computer to a time when the message is received by another 
computer. At a minimum this is the propagation delay, or the 
time it takes physically from sending to arriving in the absence 
of transmission errors and other system delays. An upper bound 
for this delay can not be given in general. It can depend on the 
amount of communication and computing going on in parallel in 
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the system, on the possibility of transmission errors (re- 
transmit tance of the message in TCP/IP) , and other random 
events, for example, page faults, process switches, and 
establishment of new communication routes) . As shown in Fig. 2, 
5 the distribution has a maximum density at a point between the 
minimum and the median delay, with a long thin tail to the 
right. In this example, 5000 roundtrip delays between two light 
H MVS processes, running on two IBM 43 81 processors connected via 
Q a channel to channel LAN, were measured (Dong, Private 
3-0 Communications June 1988) . The minimal roundtrip delay measured 
r! was min=4.22ms, and the maximal round trip delay was 

::{ 

L, max=93.17ms. 95% of all roundtrip delays measured, were smaller 
M than 5.2ms and therefore a median of 4.48ms and an average 
Q roundtrip delay of 4.91ms result. 
15 The delay has been found to be environmentally dependent. 

The delay can, under some circumstances, be approximated by a 
truncated normal distribution. Three main types of distributions 
were obsevered by Elteto and Molnar, related to rare, medium and 
busy traffic (see Figs. 3, 4 and 5) . 
20 Here again, the minimum roundtrip delay is determinated by 

propagation delay, which is in this case 64ms. The main 
difference between the distributions given in Figs. 3, 4 and 5 
is the number and location of local maxima. The busy case can be 
approximated by a truncated normal distribution, if roundtrip 
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delays beyond 600 ms are ignored. Therefore only a small number 
of "bad" values had to be deleted in this measurement (4 out of 
543) . Though this fact does not prove clearly that the 
originating distribution is a normal one, it is implied. This 
5 implication can be used to estimate the roundtrip delay, but as 
the TCP traffic differs from those traffic types, e.g. road 
traffic, arrivals at a telephone excchange, and customer 
t demands, it needs to be inferred, that the currently governing 
H distribution changes suddenly, because of the properties of the 
51.0 TCP connection. 

j?j Referring to Fig. 6, according to an embodiment of the 

U present invention at least two messages are needed to 
£ synchronize two computers, 601 and 602, connected to each other 
O via a TCP/IP socket connection. The client 601 stores at point T x 
15 its current local time. The client 601 sends the message "sync- 
request" 603, containing the time T 2 , to the server 602. As it 
takes the message an unknown amount of time t si (=transfer time) 
to be transmitted from the client 601 to the server 602, it 
arrives at time T 2 at the server 602. T 2 is the local time of the 
20 server and therefore independent from the local clock of the 

client 601. On arrival of the "sync -request" 603 at the server 
602, the server 602 stores the time T 2 . The server 602 sends a 
"sync- response" message 604 back to the client 601, which 
includes the times T lf T 2 , and T 3 , which is the time the server 
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602 retrieved from his local clock before sending the "sync- 
response" message 604 to the client 601. The transmission time 
of this message from the server 602 to the client 601 in this 
case then is t S 2- 

5 Because t S2 depends on how it is routed, this time need not 

be the same as t sl . Upon arrival of the message at the client, 
the client stores the local time in the variable T 4 . 
An offset can be estimated as follows: 

jjj_0 ©o : comparitive real offset of two machines 

Qj 0 e : estimated offset 

:., 5 : calculated roundtrip-delay 



where : 



r 2 =i;+0 o +^ 



Equation 3 



n tr Equation 4 

i 3 "" X 4 + ^0 l S2 

a - J -JJ Equation 5 

b = T 3 -T 4 Equation 6 

S-.-b-h-Td-to-T,)-!^ Equation 7 

20 ^ _ a + b _ 2G 0 +(t sl -t S2 ) ^ Equation 8 
only if 8 is near the minimum roundtrip delay (=propagation 

delay), otherwise, no estimation is possible, and therefore no 
synchronization is possible. 
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As only two equations, Equation 3 and Equation 4, result 
from the above measurements (see Fig. 6) , and these equations 
contain three unknown variables, 0 O , t S i and t S2t it may not be 
possible to calculate the value for 0 O exactly. Therefore, the 
value for 0 O can be estimated and is stored in the value 0 e (see 
Equation 8) . Let t S i=m_d+oc; (a>0) and t 52 =^_d+p/ (P>0) be the real 
time delays, whereas m_d is only the propagation delay, the 
maximum error of estimating © 0 results in: 



From Equation 9, it can be infered that, as the round trip 
delay decreases the client's error in reading the server's clock 
also decreases. Thus, if an error smaller than e shall be 
achieved, every attempt measuring a roundtrip delay larger than 
2u needs to be discarded, where: 



The closer u is to m_d, the better is the client's reading 
precesion. 

Let p be the probability that the client observes a 
roundtrip delay greater than 2u. Therefore, the probality that 
the synchronization attempt fails is p. Thus, there exists a 



maximum error: 




Equation 9 



u 



= is + m _d)~ (2e + m_d) 



Equation 10 
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trade off between the achievable precision, and the probabilty 
(\-p) of success. The better the desired precision, the smaller 
is the probabilty of success. 

This probability can be improved by allowing several 
5 successive synchronization attempts. To achieve a certain degree 
of independence between two consecutive reading at temps, they 
can be seperated by w units of time (2 seconds up to 5 seconds 
is a common value) . The probability for a successful 
0 synchronization attempt is: 

lio P =\-v k - Equation 11 

f* succes ¥ 5 

5 where: k is the maximum number of successive reading attempts 
allowed. For large values of k and a choice of w so that the 
different reading attempts can be assumed to be independent from 
I each other, Bernoulli's law states, that the average number n of 
15 reading attempts, necessary for a succesful synchronization is 



! ~5? 

S ' 



1 

n - ■ 



1- p Equation 12 

It can be said that there exist two major strategies to 
20 synchronize clocks of distributed computers. The first one could 
be called the "Aggressive Strategy" characterized by a roundtrip 
delay threshhold u close to the minimum round trip delay. It may 
need several reading attempts for a succesfull synchronization, 
but it achieves a high precesion. In this case the high risk of 
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not being able to synchronize the clocks to each other must not 
be forgotten. The second one, charcterized by a roundtrip delay 
threshold u close to the maximum roundtrip delay, needs much 
less messages for synchronization. Therefore, it bears a much 
5 lower risk of not being able to synchronize, but may achieve 
only poor precision. 

Most of the published clock synchronization algorithms 
U until 1989 assumed the existence of an upper limit max for the 
fi roundtrip delays. However, this assumption is not true. If the 
WLO roundtrip delay is smaller than the upper bound max, Lundulius 
W and Lynch showed 1984, that these methods cannot synchronize 
r. better than 



= 15 where: n is the number of clients which have to be synchronized. 

Considering two quartzes, oscillating with the frequencies 
f 2 and f 2 , both not equal to the desired frequency f 0 . Assuming 
the temperature to be roughly constant but not necessarily the 
same for the two quartzes, the difference in frequency between 
20 the two will be constant. 



(max-min)* 1 



Equation 13 



If the relative inaccuracies 



f q fo 
f0 



are A, and A 



the 



resulting offset after a measurement time t will be: 
offset = t-(A l -A 2 ) 



Equation 14 
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where A, represents the inaccuracy of the server. This clearly 
represents the equation of a straight line. 

The offset between to PCs with independent clocks increases 
almost linearly as can be seen from Fig. 7. The reason for this 
behavior can be found in the physical properties of quartz 
crystals. The frequency of each quartz depends heavily on the 
angle of the AT cut, the height of the cut and its current 
temperature. Because of technical and physical reasons, the 
manufacturing of quartzes with exactly the same frequency is 
fff.0 impossible. 

p] To compute the parameters for a line, for example, as shown 

H in Fig. 7, a linear regression can be used, which fits a 

^ straight line: 

W 

O Equation 15 

y, y = a + b-x 

15 to a given collection of N measurement pairs (y,;*,) with minimum 
mean square error. Thus, the values to be determined are a and 
b. The parameters a and b can be chosen so that the sum of the 
squares of the deviation from the straight line is minimized: 

2Q LQS^htmess^ytcompfl Equation 16 

With Equation 15 it can be concluded: 

LQS = ±Mmess)-(a + b. Xi )Y Equation 17 

1=1 
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This sum can be minimized, for example: 
dLQS = 6LQS =Q 

da db Equation 18 



Given Equation 17 and Equation 18, the values of a and b can be 
determined to be: 



N • S„, ~ S ' S 

£_ xy y — £_ Equation 19 



=3.0 



20 



S v — b ' S Y S ' S v S • S 
a=- J - - = — — y - j—^- Equation 20 



N 



N 



Equation 21 



15 ^v^X^' Equation 22 



Equation 23 



and S = Y J x r y i 

" m Equation 24 



According to an embodiment of the present invention, a 
method for estimating an offset by using linear regression is 
shown in Figs . 8 and 9 . 
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For the first thirty values, a conventional synchronization 
method, for example, one based on Flaviu's method, can be used, 
because about thirty values are at least needed to make 
reasonably precise predictions of future offsets with linear 
5 regression. To synchronize as soon as possible in order to start 
the collaborative video discussion, the method does not 
initially collect the thirty values. The initial collection of 
the values takes about 20 minutes. In this time the method can 
achieve a precision of at least 20ms. As soon as thirty offset 

140 values are measured, the parameters a and b of the regression 

ffl 

m line are computed. From that time on every 5 seconds a new 
7' estimate of the current offset is computed, by using the line, 
fU whose parameters a and jb were computed before. If the value of 
5 the current offset is smaller than 10ms, the synchronization 
**15 interval is increased by 5 seconds, and no further action is 
taken. Only if the current offset is bigger than 10ms, the 
synchronization interval is decreased by 5 seconds, and the 
estimated offset is replaced by the current offset. 

The parameters a and jb of the regression line are always 
20 computed from the 30 most recent values of the measured offset. 
Thus, again a sliding window is implemented, so that changes in 
the drift rates (because of possible temperature changes) of the 
two quartzes are taken into account. Thus, a method according to 
an embodiment of the present invention can achieve remarkable 
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improvements in the precision, as well as in the number of 
messages exchanged between the time server and the time client. 
As shown in Figs. 10 and Fig. 11, one computer, the time server, 
was connected directly to the 100Mbps LAN, and the time client 
was connected to the 100 Mbps LAN via a 19.6 kbps dial in 
connection. In the first 20 minutes, the first synchronization 
approach is used, so that the client can immediately synchronize 
to the server with a precision of at least 20ms after the 
connection is set up. When the first 3 0 offset value have been 
measured, the linear regression approach takes over. Now the 
maximum offset decreases to 5ms. As can be seen in Fig. 10, the 
synchronization interval also starts to increase continuously, 
even beyond the lOh run time. The synchronization interval has 
now reached a value of nearly 10 minutes, but the precision is 
still better than 5ms. In the first approach, the precision is 
about 10ms, and the synchronization interval did not reach 
values longer than about one minute. 

Referring to Fig. 8, a method for clock synchronization 
across a network includes waiting for a socket -connect ion 
between at least two processors across the network 801. Upon 
creation of the socket-connection the method collects statistics 
about the connection, including a measurement of roundtrip delay 
802. The method determines a roundtrip -del ay threshold 803, 
where the probability of the round-trip delay being greater than 
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the threshold is about 0.5 and the probability of the round-trip 
delay being less than the threshold is about 0.5. The method 
determines a current round-trip delay and an offset 804. The 
current round-trip delay is added to the statistics, and a new 
5 threshold is determined 805. The method determines whether the 
current roundtrip delay is greater than the new threshold 806. 
If so, the method determines whether thirty round-trip delays 
fr* have been determined 811. Upon determining thirty round-trip 
3 delays the method enters a synchronization routine. The method 
^.0 may take over from another routine. If thirty delays have not 

5 been determined the method determines a current round-trip delay 

s y 

l, and an offset 804. Upon determining that the current threshold 

L2 is not greater than the new threshold 806, the method determines 

2 

□ whether the offset is greater than an offset threshold 807. If 
: 15 the offset is greater than the offset threshold the method 

adjusts a local clock 810, decrementing by an update - interval . 

If the offset is less than the offset threshold, the method 

adjust the local clock 808, incrementing by the update - interval . 

After adjusting the clock, 810 or 808, the method determines a 
20 linear regression 809, and continues to determine a current 

round-trip delay and an offset 804. 

Referring to Fig. 9, a method according to an embodiment of 

the present invention is shown for the linear regression, 809 

Fig. 8. A current synchronization time is set 901. A flag is set 
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to TRUE 902 as a control mechanism. Upon determining that the 
flag is TRUE 903, the method determines whether the number of 
measured offsets is greater than thirty 904. Upon determining 
that the number of offsets is greater than thirty, the method 
5 removes an oldest offset from a list of offsets and adds a 

current offset to the list 905. The method determines parameters 
a and b of a regression line from the list of offsets 906, for 
example the thirty youngest offsets. Upon determining that the 
O number of measured offsets is not greater than thirty, the 
KLO method adds the current offset to the list 907. The method 
5 estimates the current offset using the regression line 908. The 
m current synchronization time is incremented and the flag is set 
£j to FALSE 909. The method determines whether the current 
m synchronization time is greater than an update -interval 910, and 
2l5 if not, inspects the flag 903, and upon determining the flag to 
be FALSE, estimates the current offset using the regression line 
908. Upon determining the current synchronization time to be 
greater than the update -interval 910, the method measures a 
current roundtrip delay and offset, 804, Fig. 8. 
2 0 According to an embodiment of the present invention, a 

method can be used for conferencing over TCP/IP networks, and 
especially for collaborative discussion of video content. The 
method takes into account that mobile clients join a video 
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discussion conference, so that the bandwidth necessary for clock 
synchronization is reduced as far as possible. 

One of the main advantages of the present invention is that 
the needed bandwidth for synchronization can be reduced, because 
5 the period between two synchronization attempts is more than 10 
times longer than in, for example, Flaviu's approach. Only every 
ten minutes 32 bytes are sent in each direction. In Flaviu's 

U approach, about every minute this synchronization message is 

Q 

0 exchanged. If the synchronization interval can be increased, the 
y.0 processing load on the time server can be reduced. 
W According to an embodiment of the present invention, a more 

L precise synchronization can be achieved as compared to the prior 
U art. A method according to an embodiment of the present 
Q invention takes into account the properties of quartz 
15 oscillators. 

Thus, a method can deal with a frequency deviation of up to 
about +/- 2-10' 4 , plus some margin for thermal effects, because 
the "master" may also have an inaccuracy of +/-10" 4 . In summary, 
the algorithm using linear regression to estimate the offset, 
20 gives a more persuasive impression. 

Having described embodiments for a method for clock- 
synchronization in distributed systems, it is noted that 
modifications and variations can be made by persons skilled in 
the art in light of the above teachings. It is therefore to be 
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understood that changes may be made in the particular 
embodiments of the invention disclosed which are within the 
scope and spirit of the invention as defined by the appended 
claims. Having thus described the invention with the details and 
particularity required by the patent laws, what is claimed and 
desired protected by Letters Patent is set forth in the appended 
claims . 
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